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Contributions to the energy levels in light muonic atoms and, in particular, to the Lamb shift 
fall into a few well-distinguished classes. The related diagrams are calculated using different ap- 
proaches. In particular, there is a specific kind of non-relativistic contributions. Here we consider 
such corrections to the Lamb shift in order a^m^. These contributions are due to free vacuum 
polarization loops as well as to various effects of light- by- light scattering. The closed loop in the 
related diagrams is an electronic one, which allows a non-relativistic consideration of the muon. 
Both kinds of contributions have been known for a while, however, the results obtained up to date 
are only partial ones. We complete a calculation of the a^m^ contributions for muonic hydrogen. 
The results are also adjusted for muonic deuterium and muonic helium ion. 
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PACS numbers: 31.30.jr, 12.20.Ds 

Recent progress of the PSI experiment on the Lamb 
shift in muonic hydrogen [1] has attracted interest to 
theory of the Lamb shift in light muonic atoms. Their 
study can provide us with information on certain nuclear 
structure effects with accuracy that is not available in 
any other experiment. 

To obtain such data, one has to be able to separate 
quantum electrodynamics (QED) effects from the nu- 
clear structure effects, and for this purpose an adequate 
QED theory providing high accuracy is required. Contri- 
butions to the energy levels in light muonic atoms and, 
in particular, to the Lamb shift fall into a few different 
well-distinguished classes. A specific theory stands be- 
hind each of them. There are corrections, the evaluation 
of which is identical for hydrogen and muonic hydrogen, 
and corrections that are specific for muonic atoms. The 
latter involve a certain part of QED, recoil effects and 
effects of the finite nuclear size. An important class of 
such specific contributions, which, in fact, also include 
the dominant term for the Lamb shift, is due to non- 
relativistic physics. 

We remind that atomic momenta in light muonic atoms 



Zam., 



l.bZrrie are compatible with the electron 



mass, while the atomic energy ~ {Za)'^m^ c^ O.OlZ'^me 
is much smaller than the electron mass. (The relativis- 
tic units in which h = c = 1 are applied throughout 
the paper.) Such an environment produces an important 
sector of corrections, which deal with a non-relativistic 
bound muon, while the QED effects are present only 



through the closed electron loops. Meanwhile the Comp- 
ton wave length of electron Ag = l/rrie determines the 
radius of the effective interaction induced by this kind of 
diagrams. The loops may be for either free- loop vacuum- 
polarization (VP) effects, related to the Uehling, Kallen- 
Sabry potential and higher-order diagrams, or the to 
light-by-light (LbL) scattering contributions. 




jyXjyXjyX 



(a) 



TVX nX 



(b) 




(c) 



* Electronic address: savely.ka rshenboim@mpq.mpg.de| 



FIG. 1: Characteristic diagrams for three basic contributions 
of light- by- light scattering effects to the Lamb shift in muonic 
hydrogen in order a^m^. Here, N stands for a nucleus, which 
may be a proton, a deuteron etc. The horizontal double line 
is for the muon propagator in the Coulomb field. 

The VP leading term is of order a{Za)'^m^ and it has 
been known for a while, while the second-order VP con- 
tribution (of order a^{Za)'^m^) was calculated with ap- 
propriate accuracy for muonic hydrogen in [2] only rela- 
tively recently. 

The accuracy of current and planned experiments re- 
quires a complete theory of non-relativistic contributions 
to the Lamb shift in order a^m^. The LbL contributions 
are depicted in Fig. [TJ while the vacuum polarization di- 
agrams are present in Fig. [2l Both kinds of contribu- 



tions have been known for a while, however, the results 
obtained up to now were only partial ones. In particu- 
lar, in the case of muonic hydrogen the contribution in 
Fig. ^ has not yet been calculated, while there are also 
some questions [3] about applicability of the so-called 
scattering approximation applied in [4] to evaluate the 
contribution in Fig. ^. 

Certain LbL contributions have specific names. The 
first one in Fig. [T] is a so-called Wichmann-Kroll contri- 
bution and it was calculated for muonic hydrogen with 
sufficient accuracy in [3|, Q. It was also reproduced in 
[4]; we also confirm this contribution. For muonic deu- 
terium and muonic helium-4 ion the results have been 
obtained in p, [3 and we confirm the deuterium result [6| 
and obtain a result for a muonic helium ion 



A^ = -0.0198(4) meV , 



(1) 



which is consistent with —0.02 meV of [7], but more ac- 
curate, and strongly disagrees with +0.135 meV of [8|. 

In our calculation we used approximations for the 
Wichmann-Kroll potential in the form 

VwK{r) = ^i^^—x 0.361662331 



X exp 



0.3728079X 



\/4.416798x2 + 11.39911X + 2.906096 



as discussed in \J\ (see also [95) and 

1. 528-0. 489x 



VwK{r) 



(Za)^10~^ 1. 374^3 + 1. 41a^2^2. 672; 



^+1 ' ^<1 



0.207x^+Q.367g;-0.413 



X > 1 



as considered in [10] and [2]. Here x = rrier. The results 
are consistent. Indeed, if higher accuracy is required, 
one can apply an exact expression [11| for Vwk(^) as a 
two-dimensional integral. 

The second term (Fig. ^) is called 'virtual-Delbriick- 
scattering contribution'. It has been calculated for 
muonic hydrogen in [J]. The calculation was based on 
jlOl , fl2j , where at first a scattering approximation was 
applied and subsequently a number of further approxi- 
mations was made. We remind that the scattering ap- 
proximation suggests that the external muon legs in the 
diagram in Fig. [1^ are on-shell (i.e. p^ = m^) and the 
muon propagator there is substituted for a free one, i.e., 
the kinematics is exactly the same if one calculates a re- 
lated Born scattering amplitude. Since atomic momenta 
Zam^ in light muonic atoms are compatible with the 
electron mass rrie^ the validity of such an approximation 
is questionable (see, e.g., the discussion in [3]). 

We, however, have proved that the scattering approx- 
imation is applicable within the uncertainty of order 
{Za)'^m^/me (in fractional units), which is at the level 
of about 1% in muonic hydrogen and deuterium and 
of about 4% in muonic helium. That is also correct 



for other simplifying approximations, which were made 
in the calculations for this contribution in light muonic 
atoms 0, [g, 0, [lO|, [l2|. A general idea of our evalua- 
tion is presented in Appendix O while the details of our 
evaluation are to be published elsewhere [? ]. 

Eventually, we conclude that the uncertainty of the 
method applied in 0, o, [3 is substantially smaller than 
the uncertainty of the related numerical evaluations for 
the Lamb shift correction in muonic hydrogen Q , muonic 
deuterium [6] and muonic helium-4 ion Q- 

The contribution in Fig. ^ has no specific name. 
Since any other LbL contribution has one ('Wichmann- 
Kroll term' and 'virtual-Delbriick-scattering contribu- 
tion'), sometimes it is referred to as a 'light-by-light con- 
tribution', which is somewhat confusing. 

This contribution has remained uncalculated for a 
while. Studying applicability of the scattering approxi- 
mation for the diagram in Fig. [1^, we have also managed 
to prove [? ] that this remaining contribution can be 
expressed in terms of the well-known Wichmann-Kroll 
term 



^% = ^^% 



(2) 



The uncertainty of this identity is of order {Za)'^m^/me 
(in fractional units), which is at the level of about 1% 
in muonic hydrogen and deuterium and of about 4% in 
muonic helium. 

By combining our results on the uncertainty of vari- 
ous approximations with the numerical results of other 
authors, we obtain the complete result for all LbL con- 
tributions of Fig. [TJ The result is listed in the summary 
table (Table III). 

With identity ([2]) proved and a possibility to obtain a 
result for the Wichmann-Kroll (AE^l) with high accu- 
racy for any light muonic atom, the uncertainty in the 
calculation of the complete LbL contribution now comes 
from the virtual-Delbriick-scattering term, which should 
determine the eventual uncertainty of the non-relativistic 
a^m^ term for muonic hydrogen, deuterium and helium- 
ion. 

Another major non-relativistic contribution in order 
a^rrij^ is due to vacuum polarization contributions. The 
VP terms of this order were studied for muonic hydro- 
gen in [1^. The diagrams are depicted in Fig. [21 which 
includes contributions of the first (Fig [2]2-(i), second 
(Fig[2^, d) and the third (FiglSl/) order of non-relativistic 
perturbation theory (NRPT). 

The most complicated terms are indeed related to the 
first line of Fig. [2l however, these diagrams were cross- 
checked due to their contributions to the anomalous mag- 
netic moment of muon jl5|, |l6| and we can rely on them. 

The contributions of the second line of Fig. [2l are spe- 
cific for muonic atoms and do not correlate directly with 
any calculation for the muon g—2. Those have been recal- 
culated completely independently of [14], as well as part 
of diagrams in the first line. 

We confirm the second-order terms of NRPT, while 
our result [17] for the third-order term (the last diagram 
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FIG. 2: Characteristic diagrams for free electron- VP contri- 
butions to the Lamb shift in muonic hydrogen in order a^m^. 



in Fig. [2]) disagrees with the one originally published in 
[l^ . After a correction [18] their result agrees with ours. 
In our calculations we used techniques developed while 
investigating second-order vacuum-polarization contribu- 
tions to the hyperfine structure of muonic hydrogen ^] . 



Term 


/\E^{2p-2s) 


/\E^{2p-2s) 


A£;""(2p-2s) 




[meV] 


[meV] 


[meV] 


1st order VP 


205.007 36 


227.634 67 


1665.7729 


2nd order VP 


1.658 85 


1.838 04 


13.2769 


3rd order VP 


0.007 52 


0.008 42(7)* 


0.074(3)* 


LbL (Fig. D) 


-0.000 71(15)* 


-0.000 73(16)* 


-0.005(10)* 


NR total 


206.673 02(15) 


229.480 40(17) 


1679.119(10) 



TABLE I: The non-relativist ic QED contributions to the 
Lamb shift /S.E{2s—2p) in light muonic atoms: hydrogen (H), 
deuterium (D), helium-4 ion (He), which includes VP con- 
tributions of the first {a{ZaYm^)^ the second (a^(Za)^m;^) 
and the third a^{ZaYm^ (Fig. [2]) order as well as a complete 
LbL contribution. The latter is a sum of contributions of order 



a{Za) 



(Fig.IIh), a'iZafm^ (Fig.©) and a^'iZaf 



(Fig. [T]:). Results marked with • are obtained in this work. 

The diagrams in Fig. [2] were also discussed in various 
papers in the context of the Lamb shift in muonic deu- 
terium [61] and muonic helium-4 ion [8]. For this purpose 
a part of the contributions was recalculated there. 

Here, we reevaluated all the vacuum-polarization con- 
tributions and the results are listed in Table HI Com- 
paring with the paper mentioned we have to acknowl- 
edge that our results are not in complete agreement with 
theirs. 

In [6| only contributions of Fig. [2]2 and h were directly 
calculated for muonic deuterium and the results for them 
agree with ours. However, that result were not applied 
there but instead the muonic-hydrogen result was 're- 
scaled'. That was achieved by assuming that the result 



for muonic hydrogen, presented in [14] in the form of 

/^E{2p - 2s) = Cs (-Y {Zafrrir , (3) 

where rrir is the muon reduced mass, can be directly 
applied for the muonic deuterium. That has not been 
claimed in jlj and is indeed incorrect and the value of 
the coefficient |T3, [Ilj 



Cs =0.118 680(12) 



(4) 



is valid only for muonic hydrogen (cf. with Cs = 
0.120 045(12) from [14], which needs a correction ^ 
as explained above). The related values for other light 
muonic atoms, obtained here. 



Cs 



0.1262(11) , 

0.270(17) , 



for /iD , 
for /iHe^ . 



(5) 



obviously differ from ^. 

The muonic-helium-4 paper [8] lacks a complete result 
and only a part of diagrams of Fig. [2] were recalculated. 
Our results are not in fair agreement, and in particular we 
strongly disagree in the contribution of Fig. [2^ for muonic 
helium. In the Wichmann-Kroll contribution (Fig.[T]2) we 
also strongly disagree with the result [8], while we agree 
with the result [7], for which we obtain higher accuracy 
(see P). 

Finally, we summarize in Table |T] a complete theory of 
non-relativistic QED contributions to the Lamb shift in 
muonic hydrogen, deuterium and helium ion up to the 
order a^m^ (see Q for references to the calculation of 
the low-order corrections). 
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Note Added: After submission of this work, we have 
calculated the LbL contribution within a static-muon ap- 
proximation (see Appendix A). The preliminary results 
for AEib, which are 0.00115(1) meV for /iH, 0.00124(1) 
meV for /iD, and 0.0114(4) meV for /iHe+, are some- 
what below the former results [4, 6, Y], but still are in 
fair agreement with them. These more accurate results 
will be reported in detail in a future publication [13]. 



Appendix A: On approximations for the LbL 
contribution 



We have proven a kind of theorem [? ] that the 
diagrams in Fig. [1] can be calculated in light muonic 



atoms (for simplicity we consider further muonic hydro- 
gen) within the static-muon approximation, in which the 
complete muon-line factor shrinks to 



-^(q) 



/ 



(2^ 



**(p)*(p + q), 



where q is the total momentum transfer to the muon line, 
^ is the wave function and the error is of the order of 

The scattering approximation [j, Isl agrees with the 
static-muon approximation within the same uncertainty. 
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FIG. 3: Reduction of diagrams with free/bound fermion prop- 
agators to static- fermion diagrams (i.e. to the static-proton 
approximation (= the external-field approximation) (left) and 
to the static-muon approximation (right). 

The proof of the 'theorem' is to be presented elsewhere 
[? ] and here we explain its main idea. The proton 
and muon are both non-relativistic (NR) particles. Their 
descriptions are very similar. 

Let us first compare the related diagrams with free 
NR fermion propagators (cf. Fig. [3]). The expressions 
for the muon and proton lines are identical. Apparently, 
one can approximate the proton line within the external 
field approach or, which is the same, by a static proton. 
To transform the complete NR expressions to the static- 
proton case, we have to neglect the proton kinetic energy 
in the proton propagators. Once that is done, after a 
chain of identities we should arrive at the external-field 
approximation. 



A reason to neglect the energy is the fact that a char- 
acteristic energy, related to a particle of the mass M, 
is Em ^ 7^/^7 where 7 = Zam^ is a characteristic 
atomic momentum. One can prove that we can expand 
using small parameters Em/^^c and Em/i- However, 
in muonic hydrogen they are of the same order since 
7 ~ TTig. The parameter Em/'^c (and the related er- 
ror) differs indeed for a muon and proton, but it is small 
for both. 

Thus as long as the muon propagator is a free one 
there is no difference in proving that we can apply the 
static-muon approximation and the static-proton approx- 
imation (see Fig. [3j). 

Meantime, in reality the situation for a muon and a 
proton is somewhat different. The muon characteristic 
momentum is of the same order as the electron mass 7 ~ 
me, and we should treat it as a bound one. 

The NR Coulomb Green function of a muon includes 



^c(^,p,pO 






|a(p))(a(pOI 



^ E-Ex^i^ ' 
a summation over all intermediate states A of continuous 
and discrete spectrum, involving energy of the intermedi- 
ates. The characteristic energy of an intermediate state is 
indeed of order of magnitude of the atomic bound energy 
Ex ~ 7^/^^/i and we can neglect it, as we already did 
in the case of the free propagators. After that, the sum 
over intermediates shrinks to the unity operator and the 
Coulomb propagator becomes equal to a free one with 
the kinetic-energy term p^/2m^ neglected 



^E 



|A){A| 



E-Ex + iO 



-> ^ 



E 



|A)(A| _ I 
E^iO E^iO 



[1] R. Pohl, F. Nez, A. Antognini, et al., Nature, (2010); to 

be published. 
[2] K. Pachucki, Phys. Rev. A 53, 2092 (1996). 
[3] M.L Fides, H. Grotch and V. A. Shelyuto, Theory of 

Light Hydrogenic Bound States, Springer Tracts Mod. 

Phys. 222 (Springer, Berlin, Heidelberg, 2007). 
[4] F. Borie, Phys. Rev. ATI, 032508 (2005). 
[5] M. L Fides, H. Grotch and V. A. Shelyuto, Phys. Rep. 

342, 63 (2001). 
[6] F. Borie, Phys. Rev. A 72, 052511 (2005). 
[7] F. Borie and G. A. Rinker, Phys. Rev. A18, 324 (1978). 
[8] A. P. Martynenko, Phys. Rev. A76, 012505 (2007). 
[9] P. Vogel, At. Data Nucl. Data Tabels 14, 599 (1974). 
[10] F. Borie and G.A. Rinker, Rev. Mod. Phys. 54, 67 (1982). 
[11] K. Blomquist, Nucl. Phys. B48, 95 (1972). 
[12] F. Borie, Nucl. Phys. A267, 485 (1976). 



[13] eprint [a?Xi7: 1005.48801 

[14] T. Kinoshita and M. Nio, Phys. Rev. Lett. 82, 3240 

(1999). 
[15] P. A. Baikov and D. J. Broadhurst, in New Computing 

Techniques in Physics Reasearch IV, edited by B. Denby 

and D. Perret-Gallix (World Scientific, Singapore, 1995), 

pp. 167-172; eprint arXiv:hep-ph/9504398 
[16] T. Kinoshita and M. Nio, Phys. Rev. 60, 053008 (1999). 
[17] V. G. Ivanov, F. Yu. Korzinin, and S. G. Karshen- 

boim, Phys. Rev. D80, 027702 (2009). See also: eprint 

arXiv:0905.4471 
[18] T. Kinoshita and M. Nio, Phys. Rev. Lett. 103, 079901 

(2009) (errata). 
[19] S. G. Karshenboim, F. Yu. Korzinin, and V. G. Ivanov, 

JFTP Letters 88, 641 (2008). 



